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Abstract. Quasi-stationary, or metastable, states play an important role in two-dimensional tur- 
bulent fluid flows where they often emerge on time-scales much shorter than the viscous time scale, 
and then dominate the dynamics for very long time intervals. In this paper we propose a dynam- 
ical systems explanation of the metastability of an explicit family of solutions, referred to as bar 
states, of the two-dimensional incompressible Navier-Stokes equation on the torus. These states 
are physically relevant because they are associated with certain maximum entropy solutions of the 
Euler equations, and they have been observed as one type of metastable state in numerical studies 
of two-dimensional turbulence. For small viscosity (high Reynolds number), these states are quasi- 
stationary in the sense that they decay on the slow, viscous timescale. Linearization about these 
states leads to a time-dependent operator. We show that if we approximate this operator by drop- 
ping a higher-order, non-local term, it produces a decay rate much faster than the viscous decay 
rate. We also provide numerical evidence that the same result holds for the full linear operator, 
and that our theoretical results give the optimal decay rate in this setting. 



1. Introduction 

It has been observed, both numerically and experimentally, that solutions to the nearly-inviscid 
two-dimensional incompressible Navier-Stokes equation on the torus will rapidly approach certain 
long-lived quasi-stationary states [H [131 02] • Evidence also suggests that these quasi-stationary, 
or metastable, states are connected with stationary solutions of the inviscid Euler equations. In 
|18j . the authors argue that certain maximum entropy solutions, determined using numerical and 
analytical techniques, of the inviscid Euler equation are the most probable quasi-stationary states 
that one would observe. They refer to these states as dipole and bar states and confirm their 
predictions by numerically analyzing the Navier-Stokes equation and demonstrating that most 
solutions do indeed converge rapidly to one of these two classes of states. Furthermore, in [3], the 
authors consider a stochastically forced Navier-Stokes equation, and observe both the dipole and 
the bar states as statistical equilibria. 

There are explicit exact solutions of the Navier-Stokes equation that qualitatively match the bar 



and dipole states. (See equations (1.2) and (5.1), below.) The bar states are often referred to 



as Kolmogorov flow |2J, and the dipole as a Taylor-Green vortex [121 Pg- 172], but we adopt the 
terminology of |18j when referring to these exact solutions. In this work, we provide analytical and 
numerical results indicating that there is a class of initial data for which solutions to the linearized 
equation rapidly converge to the bar states. More precisely, for small viscosity < v <C 1, we prove 
that for a natural approximation to the linearized operator, there exists an invariant subspace for 



Date: May 18, 2012. 



2 



MARGARET BECK AND C. EUGENE WAYNE 



the linearized equation in which solutions converge to the bar states at a rate that is 0{e~^ 1 ), 
which is much faster than the decay of the bar states themselves, which is 0{e~ ut ). Furthermore, 
we present the results of numerical computations that indicate that this decay rate is optimal, not 
just for our approximation but also for the complete linearization of the Navier-Stokes equation 
about the bar states. We do not analyze the dipole solutions here - see ^5] for more information 
about these states. 

This work can be viewed as an extension of our prior research on metastable behavior in Burgers 
equation, [lj. In that work we showed that corresponding to each long-time asymptotic state of 
Burgers equation there is a one-dimensional, invariant manifold of viscous iV-waves which attracts 
nearby solutions at a rate much faster than the viscous time scale and which governs the behavior 
of solutions for very long times, before they reach their ultimate, asymptotic state. Furthermore, 
this manifold is tangent, at the fixed point representing the asymptotic state, to the eigenspace of 
the eigenvalue of smallest real part. For the two-dimensional Navier-Stokes equation on the torus, 
the asymptotic state is the rest solution in which the velocity of the fluid is zero. Linearizing about 
this state, the eigendirections corresponding to the smallest non-zero eigenvalues correspond to the 
bar states, and these invariant manifolds (in this case, actually invariant subspaces) form attractive, 
invariant manifolds of metastable states, just as in the case of Burgers equation. 

It is interesting to note, however, that the mechanism causing the separation in times scales for 
Burgers equation and the bar states of the Navier-Stokes equation appears to be different. Lin- 
earization about the bar states produces small eigenvalues that scale with that square root of the 
viscosity parameter. In contrast, linearization about the asymptotic states for Burgers equation 
leads to an operator whose eigenvalues are independent of viscosity, and the separation in time 
scales appears to result from exponentially large coefficients in the eigenfunction expansion. Both 
phenomena, however, are connected with the fact that neither linearized operator is self-adjoint. 

Consider the two-dimensional incompressible Navier-Stokes equation on the torus, 

dt.u = vAu — u • Vu + Vp, V • u = 0, 

where — tt < x < tt, —tt < y < n, u : T 2 — > M?, u(— n, y, t) = u(ir, y, t) and u(rc, — tt, t) = u(x, tt, t) 
for all t > 0. The pressure p(x, y, t) is a scalar valued function that also satisfies the periodic 
boundary conditions, and the viscosity v is assumed to be small: < v <C 1. The fluid vorticity oj 
is defined in terms of the velocity u via ui = V x u = d x U2 — d y u\. Taking the curl of the above 
equation, we find the two-dimensional vorticity equation 

(1.1) dtOJ = vAuj — u • Vw, 

where the scalar-valued function ui also satisfies periodic boundary conditions. We assume without 



loss of generality that solutions to this equation have zero mean - equation (1.1) conserves mass 
and, in fact, the periodic boundary conditions for u imply that ui has zero mean. In this case, one 
can determine the velocity from the vorticity using the Biot-Savart law as follows. Given a periodic 
function oj £ L 2 (T 2 ), we can write it as its Fourier series 



u(x,y) = 2^ Qj{k,l)e l 
k,iez 
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where 

oj(k,l) = ~ [ Lo(x,y)e-^ kx+ly Uxdy. 

4VT J [—TTjTr] X [— 7T,7r] 

Using the fact that w = Vxu and the incompressibility condition V-u = 0, we have the Biot-Savart 
law 

= ^Tpih-Wfal), (k,l) (0,0). 
Alternatively, one can obtain this relationship using the fact that, if tp is the stream function, then 
Atp = w and u = V ip = {—d y tp : d x ip). Since w has mean zero, 

u= (-d y A~ 1 w,d x A- 1 w). 



The explicit family of solutions of (1.1) that qualitatively matches the maximum entropy solutions 



found in [181 Figure 2c] and referred to as "bar states" has vorticity given by 
(1.2) uj b ' m (x, y, t) = e~ vmH cos(mx), u b ' m (x, y, t) = e~ vmH sin(mx) 

and the associated velocity solutions are 

.6,m/„ „, j.\ A um 2 t I ^ l ~h,mf ^ — / 





u°' m (x,y,t) = -e~ vm f . , u°' m (x,y : t) 

smmx I m \cosmx 

Note that the bar states can also be defined to be functions of y, rather than x, and any linear 
combination of x-bar states, or any linear combination of the y-bar states, is again an exact solution. 

We will focus primarily on these states when m = 1 and denote ui = u^' 1 . When < v -C 1, 
these states are quasi-stationary, because they evolve in time only on the very long timescale, 
0{e~ vt ). They are stationary states for the inviscid Euler equation (i.e. when v = 0). Due to 



incompressibility and periodic boundary conditions, for any solution of (1.1) we have 

where we have used Poincare's inequality to obtain the final inequality. Thus, the slow decay rate 
of the bar states occurs on what is arguably the natural timescale for the system, determined by 
the viscosity v. 



Consider a solution of ( 1.1 ) of the form ui = aw + w, where the constant a represents the amplitude 



of the bar state around which we will linearize. Since equation ( 1.1 ) conserves mass, we will assume 
that uj, and hence also w, has mean zero. The perturbation then satisfies 

(1.3) d t w = C{t)w - u • Vu> 

C(t)w = vAw - ae~ ut [sinxd y (l + A" 1 )]?/;, 

where u is the velocity corresponding to w via the Biot-Savart law. 



Our goal is to show that there exists an invariant subspace for the linear part of equation (1.3) 
on which solutions decay very rapidly, at a rate that is 0(e~^ t ). See Theorem [l] There are a 
couple of issues which have to be overcome to obtain this result. First, ignoring for the moment 



the time-dependence of the operator in (1.3), it must be shown that the addition of the term 
sin xd y (l + A -1 ) to the Laplacian produces increased decay. In the current paper we prove this 
for an approximation to the linear operator which ignores the nonlocal part of this operator. This 
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is done using the so-called hypocoercivity properties of the operator as developed in [17]. We also 
discuss the changes that will be needed to extend this result to the full linearized operator. We note 
that a similar approach has recently been used to study the linearization of the two-dimensional 
Navier-Stokes equation around an Oseen vortex [5j. 

Second, it's clear that we can't expect decay at the rate 0{e~^' t ) on the entire space. This 

,2/ 



is because there are explicit solutions of the linear part of equation (1.3) given by e~ um i=timx 
and e~ utziziy . Therefore, we need to find a way to project off of these anomalous modes. If the 
operator was time-independent, this could be done using a spectral projection via the adjoint 
eigenfunctions associated with the eigenvalues — v and —m 2 v. However, that will not work in 
this case. Nevertheless, we are able to find an invariant subspace in which solutions are rapidly 
decaying. We first examine this subspace for the full linear operator, and then show that for our 
approximation, in which we omit the nonlocal term, the rapidly decaying modes have a simple 
expression in terms of the Fourier coefficients of the solution. 

A further consequence of the time-dependence of the operator is that we cannot expect a result 
that holds for all time. It is the term ae~ ut \sm.xd y (l + A^ 1 )] that must be producing the rapid 
decay, and this term is 0(1) only for times up to and including 0{\/v). Hence, Theorem [T] is only 
valid on a time interval of this size. 

The reason that we focus on the bar states, rather than the dipoles, is primarily that these states 
seem easier to analyze from a mathematical point of view. Although in [TH] it was found that a 
large class of initial conditions have the dipole as their quasi-stationary state, we note that there 
are certain situations in which bar states seem to be more physically relevant. For example, in [3] 
the authors find that the bar states become dominant if the torus is elongated slightly in either the 
x or the y direction. 

There has been much previous work related to the multiple time-scales in the evolution of solutions 
to the two-dimensional incompressible Navier-Stokes equation. Those most relevant for the present 
analysis have been referred to above. We do not attempt to provide a comprehensive list, but we 
mention a few other works that are in the same spirit, in the sense that they approach the problem 
from a dynamical systems point of view. In particular, in [8j [10] the authors analyze the global 
stability of the single Oseen vortex for the 2D Navier-Stokes equation on the plane using invariant 
manifolds and a Lyapunov functional, and in [7] the author provides rigorous asymptotic expansions 
for the slow motion of localized vortices, again for the 2D Navier-Stokes equation on the plane. The 
case of two identical vortices is considered in [9] . The paper of Lan and Li [JT] considers eigenvalues 
of the linearization of the Navier-Stokes equation for small values of the viscosity but does not 
examine the viscosity dependence in the way we do here, nor do they consider the relationship of 
this problem to metastable behavior. We also note that there is potentially a connection between 
the \fv scaling, found here, and Taylor dispersion [16], although we have not explored this in detail. 

The rest of the paper is organized as follows. In ^2] we identify a small set of anomalous modes that 
decay on the viscous timescale 0(e~ ut ) and define an invariant subspace containing none of these 
modes. In £j3j we prove that for a slightly simplified approximation to this linear operator, solutions 
of the associated evolution equation decay at the rate of at least 0{e~^ 1 ). We also discuss the 
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possible extension of these results to the full linear equation. Next, in ^4] we present numerical 
evidence that the decay rate derived in the prior section is in fact optimal by studying the way in 
which the eigenvalues of the linearized operator depend on the viscosity, v. Finally, we conclude 
with a discussion and mention directions for future work in §[5j 



2. Anomalous modes 



Our main goal of this paper is to show that the bar states attract nearby solutions at a rate much 
faster than expected from the viscous time scale. However, there are some "anomalous" modes 
which do decay only at the viscous rate. In the present section we identify the set of states with 
this anomalous behavior and show that it is "small" in the sense that if we expand solutions in a 
Fourier series with respect to the basis {e ifcx e J ^}, then no modes with |/| > 1 behave anomalously, 
and only a very special set of those with \l\ = 1 do. 



The linear part of equation (1.3) is given by 
(2.1) 



wt = C(t)w, C{t) 
In Fourier space, this operator becomes 



i/Aw - ae"^[sinx<9 y (l + Ar l )]w. 



£(t)w(k,l) 



-vik 2 + l 2 )w(k,l) 



-ae 



-vt 



(k-l) 2 + p 



w(k - 1,1) 



(2.2) 



{k + l) 2 + l 2 



w(k + 1,/) 



We wish to show that solutions to equation (2.1) decay at a rate which is faster than 0(e ). As 



mentioned above, there is one key obstacle preventing this from being true, in general, and that 
is the presence solutions of the form e - um2t e ± imx anc | e~ ut e ±iy , which we refer to as anomalous 
modes. Therefore, any solution whose Fourier series has d)(±m, 0)(i) 7^ and/or w(0,±l)(i) 7^ 
cannot decay faster than 0(e~ um2t ). 

In some sense these solutions correspond to the eigenvalues —vm 2 for the linear operator. How- 
ever, since the linear operator is time-dependent, these are not eigenvalues in the usual sense. 
Nevertheless, we will see they have a similar effect, and we will slightly abuse this terminology. If 
the operator were time- independent, we could simply determine the corresponding adjoint eigen- 
functions and use them to construct spectral projections and an invariant subspace without any 
anomalous modes. However, in this time-dependent case we must do something slightly different. 



Lemma 2.1. The functions e um2t e ± imx } m S N, and e ut e ±iy are solutions of (2.1). In addition, 



' are solutions of the corresponding adjoint equation dtw = C*(t)w, where C* is the L 
adjoint of C, but e~ vt e ±iy are not. 

Proof. This follows from a direct calculation, using the fact that 

C*(t) = vA + ae- ut [(l + A" 1 )^ sinx]. 
The key observation is that d y e ±mix = (1 + A~ l )e ±iy = 0. 1 
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Because the functions e - um2t e ± imx a i so solve the adjoint problem, the subspace 



■Mrapid,* : = G L2 ( j2 ) : 0) = 0, m £ 7j} 



is invariant for the linear PDE (2.1). Therefore, we can avoid these anomalous modes by working 
only in this subspace. Unfortunately, it is less straightforward to remove the other anomalous 
modes, because the corresponding adjoint "eigenfunctions" are not time-independent, as the fol- 
lowing calculation shows. 

Suppose we look for solutions to — vu) = C*u, in attempt to determine a projection off the modes 
e ±iy . In Fourier space, this equation is 



= i/[l- (k 2 + l 2 )]cb(k,l) 



I 



-ae 



-vt 



1 



1 



[w(k — 1,1) — w(k + 1, 1)]. 



k 2 + P 

As mentioned above, e ±imx are solutions. Are there others? When (k,l) = (0,±1), the above 
equation is always satisfied. Thus, without loss of generality we can take w(0, ±1) = 1. However, 
we then find for (k, ±1) = (1, ±1) 



= -vu(l, ±1) ± -ae- ut [w(0, ±1) - w(2, ±1)], 



and for (fc,±l) = (-1,±1) 

= 



±1) ± ^ae- ut [w(-2, ±1) - w(0, ±1)], 



and so on. Thus, the effect of the (0, ±1) mode will propagate out requiring nonzero, and time 
dependent, entries in other components of the adjoint eigenfunction. Formal calculations with this 
equation lead one to the following result. 

Define 



(2.3) pf := w(2j, ±1) + w(-2j, ±1), qf := w(2j + 1, ±1) - w(-(2j + 1), ±1), 



for j € N and 
(2.4) 



■Mrapid,!/ :={w£L 2 (T 2 ):pf =qf =0 for all j £ N}. 



Proposition 2.2. A solution w(t) G L 2 (T 2 ) of (g.J| ) satisfies w(0, ±l)(t) = for all t > if and 
only if w(0) £ A^ rap id,y, defined in (2.4)- 



Proof. Using equation (2.1), we find 



-v{\ 3 2 + l)pf =F 



1 



-v((2j + l) 2 + l)qf T -e^ 



(2j - I) 2 + 1 
1 



1 



3-1 



1 



(2j + 1) 2 + 1 

1 

(2j + 2) 2 + l 



(2j) 2 + l, 

As a result, the space M ra , p id,y is invariant, so if w(0) £ A4 rap id )3/ then w(0, ±l)(t) = for all t > 0. 

To prove the other direction, suppose u)(0, ±l)(t) = for all t > 0. Using the Fourier representation 
of £(i), given in (2.2), we see that the dynamics of w(0, ±1) are determined by 



3tt&(0, ±1) = -uw(0, ±1) + -Ae- Ut [w{-1, ±1) - tu(l, ±1)]. 
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Thus, in order to have w(0,±l)(t) = ^Po = for all t > 0, it must be the case that w(—l, ±1) 
w(l,±l) = —q^ = for all t > 0. Looking now at the above equations for p^ and </• , we have 



4 ± 



which then implies we must have pf = for all time. One can continue in this manner to 
show that qp = pf = 0. Another way to think about this is that, if one were to define = 
(Po j 1o > pf > it > ■■•)) then -u^ = ^(t)^, where is a tridiagonal matrix. The only solution to 
this linear equation with u^(t) = (0, 0,^ (t), (t), . . . ) is the solution that is identically zero. This 
proves the result. ■ 

As a result, we define 

(2.5) M := A^ rap id,x n M r£kpid ,y, 

which is the subspace one which one can hope to have rapid decay of solutions. In [18] the authors 
found rapid convergence to the bar states only for specific types of initial data. (More general initial 
data converge rapidly to the dipole solutions - see ^5]) Examples of such initial conditions that 
they found do indeed seem to match the condition imposed by the subspace M. . See [181 Figures 8 
&; 13]. Note, however, that in spite of the restrictions the anomalous modes impose, the bar states 
are still physically relevant due to the fact that they can be viewed as maximum entropy solutions 
of the Euler equation, as well as the fact that they seem to dominate the flow when the aspect ratio 
of the torus is changed slightly [3]. 



3. Rapid decay via hypocoercivity 



In the present section, we wish to study the decay generated by the linearized operator (2.1). In 
fact, we will study an approximation to that operator obtained by dropping the non-local operator 
(sinx)3j / A _1 , and studying the approximate equation 

d t w = C a w 

(3.1) C a w = vAw - ae~ ut (sin x)dyW . 

Note that the term we dropped is a relatively compact perturbation of C a and hence we believe 
that C a should capture the most important aspects of the linearized evolution. Furthermore we 
present numerical results on the full linearized operator in Section [4] that support this point of view. 
We discuss at the end of this section the extension of our current results to the full linear operator. 
We also note that a similar approach has recently been utilized to study the linearization of the 
two-dimensional Navier-Stokes equation around the Oseen vortex [5]. 

To state our results, note that if we wr ite the solution w = J2 k i w(k, l) e i{kx+ ^ of dO> ELS cl Fourier 



series, the operator L a does not mix different values of the Fourier index £. Thus, if we define 
(3.2) w(x, y) = J2(Mx)V ey = E ( E ^ l ) e ) ^ 
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we can rewrite 



(3.3) C a w = J2(£,2wi)e Uv 

t 

where 

(3.4) C%wt = u(d 2 x - £ 2 )w e - ia£e~ vt (sin x)w t . 

Motivated by Villani's work on hypocoercivity, |17j . we define operators, 

(3.5) B e w e = -ia£e~ ut (sin x)w e 

C e we = [d x , B l ]wi = -ia£e~ ut (cos x)we ■ 

As we will see, the fact that [d x , B ] ^ is responsible for the faster than expected decay of solutions 
of our linear equation. 

Remark 3.1. The discussion of anomalous modes from the previous section is simplified for C a . 
The operators B e and C l both annihilate any functions wq, while there are no nonzero elements of 
the kernel of B e or C for any £ ^ 0. Thus, we will show that the only anomalous modes are those 



with £ = and we have accelerated decay for all solutions of (3.1) of the form 
(3.6) w(x,y) = J2Mx)e iiy ■ 

Given the preceding remark, we define the Banach space 
(3.7) 



Ml 2 + \lui\\ d xwe\\ 2 + ^lo/j C^il 



\w\\ x < oo 



X = { w e L 2 (T 2 ) : w = 0, ^ 

Note that the I = modes are absent from functions in X. In this definition, and throughout 
the remainder of the paper, || • || will denote the norm on L^T 1 ). We will use (•,•) to denote the 
innerproduct on L 2 (T 1 ). 

Theorem 1. Given any constant r and T G [0, t/u], there exist constants K and M that are 0(1) 
with respect to v such that the following holds. If v is sufficiently small, then the solution to the 



linear equation (3.1) with initial condition w° G X satisfies the estimate 

\\w(t)\\ 2 x < Ke-^WwYx 

for all t G [0,T]. 



Hence, at least in our approximation, the bar states are more stable than one would expect from 
the diffusive decay of 0(e' vt ). In fact, if we take T = t/u, we find e~ M ^ T = e~ Mr < e~ uT = e~ T . 



(See Remark 3.5 below, for a discussion of how the constants M and K are related to T.) Therefore, 
even though this result is only valid for finite times, it still implies that, for initial data in the above 
mentioned subspace, there will be an initial period of rapid decay to the bar states, which then 
themselves decay on the diffusive timescale. 

As we mentioned in the introduction, this result is reminiscent of our prior work on metastable 
states in Burgers equation, [lj. As in that work, the bar states form an invariant manifold (in 
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fact, an invariant subspace in this case) which connects to the long-time asymptotic state of the 
system (the zero state in this case), and at that point, it is tangent to the eigenspace of one of 
the eigenvalues with smallest non-zero value. Furthermore, the preceding theorem shows that (at 
least in our approximation), nearby solutions are attracted to the invariant manifold at a rate 
much faster than the motion along the manifold, justifying the identification of these solutions as 
metastable states of the system. 

The remainder of this section is devoted to the proof of this Theorem. 



Consider equation (3.1). As we said above, we wish to use Villani's hypercoercivity method, as 
developed by [6|, and as a consequence, the properties of the operators B l and C introduced in 
the introduction of this section will play an important role. Thus, we begin, with a simple lemma 
detailing some of the properties of these operators. 

Lemma 3.2. The operators B £ and C e are bounded, anti- symmetric operators on L 2 (T l ) for any 
£ ^ 0, with bound 

(3.8) ll-B^II < \a£\e~ ut \\w e \\; \\C%\\ < \a£\ e - ut \\w £ \\ . 

Furthermore, B^C e = C^B f ~ - i.e. these operators commute. 

Proof. The proof is a simple computation using the definition of £? and C . For instance 
\\B e w e f = (a£) 2 e- 2ut J sin 2 x\w e (x)\ 2 dx < {a£) 2 e~ 2vt \\w e \\ 2 

while 

(we, (B e ue)) = a£e~ ut J w^{x){—iuf {x))dx = —a£e~ ut j (—iwe(x))ue(x)dx = —(B i we,ue) 

with a similar computation for ||(7^w)£||. ■ 
Because we wish to follow the scheme developed by Gallagher, Gallay and Nier, we we rescale the 



time by t — > vt and consider the following problem equivalent to (3.4), 



(3.9) d t wt = (<9 2 - £ 2 )w e + ^B e w e , £ ^ . 



Following [T7] and [BJ, we will measure the changes in solutions of (3.9) in terms of the functionals 
(3.10) ¥(t) = \\w e \\ 2 + a\\d x w e \\ 2 -2pRe(d x w e ,C i w e ) +j(C l w £ ,C e w e ). 



We will establish decay for solutions of (3.9) for each £ by choosing the constants a, (3 and 7 



appropriately, and then combine these estimates to prove Theorem [TJ 
Remark 3.3. Note that 

\2(3Re(d x w e ,C e w e )\ < " \\(w t ) x \\ 2 + ^ \\C £ w e \\ 2 . 

2 a 

Thus, if 

(3-11) /? 2 < ^, 

we have 

\\M 2 + f MmU 2 + l\\c*w e \\ 2 < &(t) < \\w e \\ 2 + ^IK^II 2 + ^ll^ll 2 - 
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Consider the time change of this functional. Also, for the moment, we will be working with equation 



(3.9) for a fixed value of £, so we suppress the subscript £ on the functions wg to avoid overburdening 
the notation. 

^$*(*) = ((w t , w) + (w, w t )) + a ((d x w t , d x w) + (d x w, d x w t )) 
(3.12) -2pte({d x w u C l w) + {d x w,C l w t j) +l[(C e w t ,C e w) + (C e w,C e w t )^ 

-2/3 Re(d x w, ^w) + 7 ((^w, C'w) + (C*w, ^w) 



We bound in turn each of the pairs of terms in (3.12) 



3.1. The term ((wt,w) + (w, Wt))- Using ( |3.9[ ) we can rewrite this as: 

(3.13) {{{-£ 2 + d 2 + ^B e )w,w) + (w, (-£ 2 + d 2 + J#V)J = -2£ 2 \\w\\ 2 - 2\\w x \\ 2 , 

where we used the anti-symmetry of B to eliminate terms like (w, B e w). 

3.2. The a term ((w x t,w x ) + (w x ,w x t)). Again, we rewrite this as: 
(d x {-£ 2 + d 2 x + ^B e )w,w x ) + (w x ,d x (-£ 2 + 8 2 x + = -2£ 2 \\w x \\ 2 - 2\\w xx \\ 2 

+ ^[(d x (B e w),w x ) + (w x ,d x (B £ w))} 
Using the anti-symmetry of B £ , we can rewrite the last two terms as 

[(d x {B l w),w x ) + (w x ,d x (B e w))] = (B l w x ,w x ) + ([d x , B £ ]w, w x ) + (w x ,B £ w x ) + (w x , [d x ,B l ]w) 

= 2Re(w x ,C e w) . 

Since |2 He(w x , C^w)] < 2\\w x \\ \\C w\\, we finally obtain 
(3-14) 

a \{-£ 2 + d 2 x + -B e )w,w x ) + (w x ,d x (-£ 2 + d 2 + -B e )w] < -2a£ 2 \\w x \\ 2 -2a\\w xx \\ 2 + — \\w x \\ \\C e w\\. 



3.3. The term {(d x w t ,C l w) + (d x w, C^wt)). Using ( |3.9[ ) we can rewrite this as: 

(d x w t , C e w) + (d x w, C e w t ) = {d x (-£ 2 + d 2 x + Ib £ )w, C e w) + (d x w, C\-£ 2 + d 2 x + ^B l )w) 

= -2£ 2 Re(d x w, C l w) + [{w xxx , C l w) + (w x , C e w xx )} 
+ ^[(d x (B e w),C e w) + (w x ,C £ (B £ w))] 

By Cauchy-Schwartz inequality, we can bound \(w x ,C e w)\ < \\w x \\ \\C e w\\ . Integrating by parts, 
and using the fact that [d x , C £ ] = — £? , we can write 

(w xxx , C e w) + (w x , C e w xx ) = -(w xx , d x (C l w)) + (w x , C l w xx ) 

= -2Re(w xx ,C e w x ) + (w xx ,B e w) 

= -2Re(w xx , C l w x ) - (w Xl d x B l w) 

= -2Re(w xx ,C e w x ) - (w x ,B e w x ) - (w x ,C £ w) . 
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Because B^ is anti-symmetric, the middle term has zero real part, and so it will not contribute 
anything to the time derivative of <E> (t), as we only consider the real part of this term. We can 
therefore bound this expression by 

\{w xxx ,C e w) + (w x ,C e w xx )\ < 2||w xx ||||C^u; x || + ||w x ||||C^|| 

Next, commuting d x and B e , we have 

(d x (B e w), C e w) + (w x , C\B £ w)) = (B e w x , C e w) + (C e w, C l w) + (w x , B £ C e w) + (w x , [C e , B e ]w) 

= \\C e w\\ 2 + (w x , [C e , B e ]w) = \\C e w\\ 2 , 

where the penultimate equality used the anti-symmetry of B l and the last equality used the fact 
that [B e , C e ] = 0. 

Thus, combining the above estimates, we obtain 

(3.15) -2pRe[(d x wt,C e w) + {d x w,C e w t )} < (4^+2)/3||^ x ||||C^||+4/3||u; ra ||||C^||-^||C^|| 2 . 

3.4. The 7 term ((C e w t ,C e w) + (C e w,C e w t )). Rewrite this as 

{C\-l 2 + d 2 x + ^B e )w, C e w) + (C e w, C\-t 2 + d 2 x + ^B e )w) 

= -2£ 2 \\C e w\\ 2 + {{C e w xx , C e w) + (C e w, C e w xx )^j + ^ ((C e {B e w), C l w) + (C E w, C e (B E w))^ 

We have 

(C e w xx ,C e w) + (C E w,C e w xx ) = -[(C e w x ,d x (C E w)) + (d x (C e w),C e w x )} + [(B e w x ,C e w) + {C e w,B e w x )] 

= - [{C e w x , C l d x w - B l w) + (C e d x w - B e w, C e w x )} 

+ [(B e w x ,C e w) + (C e w,B £ w x )] 
= -2{C e w x , C e w x ) + (d x C e w + B e w, B e w) + (B e w, d x C e w + B e w) 

+ [(B e w x ,C e w) + (C e w,B £ w x )\ 
= -2{C e w x , C e w x ) + 2(B e w, B e w) - [(C e w, d x B l w) + (d x B e w, C e w)} 

^(B e w x ,C e w) + (C e w,B e w x )} 
= -2{C e w x , C e w x ) - 2(C e w, C e w) + 2(B e w, B l w). 

The term proportional to u^ 1 can be computed using the fact that 

(C e (B e w), C l w) + (C e w, C e {B e w)) = {B e (C e w), C e w) + (C e w, B e (C e w)) = 0, 

where the first step used the fact that C and B commute, and the second step used the anti- 
symmetry of B . 

Hence, we have 

(3.16) j[(C e w t , C e w) + (C e w, C e w t )} < -{2l 2 + 2) 7 ||C^|| 2 - 2 7 ||C^ :r || 2 + 2 7 ||£^|| 2 . 
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3.5. The terms containing Noting that = —vC 1 , (from ( |3.5[ )), we have 



(3.17) 
while 
(3.18) 



Ant 

\2/3Re(d x w, —w)\ < 2^11^1111^1 



3.6. Estimate on the time evolution of <1> . Combining equations (3.13)-(3.16), we obtain 



j&(t) < -2£ 2 \\w\\ 2 - [2 + 2a£ 2 ]\\w x \\ 2 - 2a\\w xx \\ 2 + + 2(3(2£ 2 + 1 + ^ )]|KI|||C^[| 

+4/311^111167^11 - \{2i 2 + 2) 7 + ^ - 2 7 ^]||C^|| 2 - 2 7 ||C £ ^|| 2 + 2 7 ||i?^|| 2 . 
We now use the facts that 



4f3\\w xx \\\\C e w x \\ < a\\w xx \\ 2 + ^\\C e w x \ 

a 



2a, 



a 



<,t.l|2 



< \\w x \\ 2 + ^\\C l w 



2/3(2£ 2 + l)|K||||C^|| < l a (2l 2 + 1 + ^)||^ 2 || + + 1 + v) \\C e wf 

2 a 

to obtain 

j-$'(t) < -2^ 2 ||u>|| 2 + 2j\\B e w\\ 2 - [2 + 2a^ 2 - 1 - ^a(2£ 2 + 1 + u)]\\w x \\ 2 - (2a - a)\\w xx \\ 2 

(JUL — 



4 /3 2 i„^. „ 2 ,2^^,. « 2 2/3 2 (2^ 2 + l + z,) ... - t| . a 



'[27 " — ]||C^|| 2 - + (2l 2 + 2) 7 

a v v 



a 



2jv]\\C e w\ 



Now set a = ao\As P = Po, and 7 = 70/-^) where ao, A), and 70 will be assumed to be independent 
of v. to find 



d 



*«(t) < -2^wf + 2^||B^|| 2 -[(l-z, 



) + ao\/^^ 2 ]||^|| 2 - aoV^ll^a 



2( 7o a - 2/3 2 



+2 



||c ^ ||2 _ [ (2/3o-a 2 L ) , np2 



+ 2t 



(7o«o - 2/3 2 ) 
a v^ 



(7o«o - /? 2 ) 2/3 2 



ao 



^-27ov^]||C^ 



£„-mi2 



Hence, as long as ao, A), and 70 are chosen so that (3.11), and 

(3.19) fa > 4a 2 
hold, and if v is sufficiently small, we have 

(3.20) -^^>^C*) < — 2^ ll^ll 2 — ^ ||^^ II s — ^o^^M^^-^ II s — ||<^^^_ || 2 — MC^^^II^. 
at \V 4 Jv _'/' 



Our goal is to show that ($)'(£) < ^=<1> (i). All terms in (3.20) are of the correct form except for 



— 2^ 2 ||tt;|| 2 + ^= ||i?^i()|| 2 . We handle these terms with the aid of a proposition modeled after the 
analogous bound of [U Proof of Prop 4.1]. 
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Proposition 3.4. Ifr, T G [0,t/v], a, and \£\ > 1 are fixed, there exists a constant Mq such that, 
for v sufficiently small and for all t G [0, T], 



(3.21) 

for any u G i? 1 (T 1 ). 



1 



u. 



,2 , ftii^i^W, 



+ 



2zV 



Remark 3.5. The constant Mq can be chosen uniformly for t G [0, T], but its value will depend on 
T. The proof below indicates that Mq ~ 0{e~ uT ). This makes sense: since \\C l \\ ^ ast — > oo, 
we must have Mq — > as t — > oo, as well. This implies we may take T = 0(1/ v), since then 
Mq = 0(1), in the statement of this proposition, and therefore also in the statement of Theorem^ 



Before proving this proposition, let's look at its consequences for equation (3.20). Since ||-E> w|| < 
ale~ vt \w\ < al\\w\\, equation (3.21) implies 



(3.22) 

We now define 
(3.23) 



v 



1 



^T.l|2 



2z/ 



\C'w 



<-( 



2/2 



\W\ 



~Mf 



which implies that the right hand side of (3.22) is less than or equal to 

2 7o a 2 * 2 



w 



Next, to insure that (3.19) holds, we set 
(3.24) 



2 _ 0o _ 4 7o 2 a 4 l 2 



Finally, we must check that (3.11) is satisfied. Given the definitions of ctQ and (5q, this inequality 
will hold provided 



7o > 



4/3o 2 _ c g3/2_^167o 2 « 4 ^- :2 



CtQ 



Ml 



Therefore, combining this with equation (3.23) and (3.24), we take 

( 3 - 25 ) 7o = — t,;,,,,,, , <">,, 



M, 



1/2 



' 64V2a 3 |£| 3 / 2 ' 
With these definitions, we obtain: 



32 v / 2a|£| 1 /2 



Bo 



Mq 



32(16)a 2 |£|' 



(3.26) < -^ 2 \\wf 



2-fQaH 



2/2 



1 



\W\ 



Wx\\ 2 ~ CtQy/v\\w x 



^IIC^II 2 -^ 



\C e w x \ 



where ctQ, Po, and 70 are defined in (3.25). 

We wish to show that (3>^)' < —(M/^/v)<& for some M, which will hold if <£> / \fv is less than or 
equal to —1/M times the right hand side of (3.26). Remark 3.3 implies that 
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which needs to be less than or equal to 



1 

M 



2/2 



2t\\w\\ + 



2 7o a^ 



1 



A), 



Mr + olhM + «ov^Na«ll + —lie ^11 + 

o f 



70 icVll 2 



Using equation (3.25), we find this will hold if we chose M such that 

.r 3 / 2 

M < 2-f a 2 £ 2 



Mn 



M < 



M < 



32^/2a 
1 8V^av1^[ 



12a 3VMo 

2/Jo = 24a 

3 7o ' 



proof of Proposition 



3.4 



Since \£\ > 1 on the sub space M, we can take M = 0(Mq' ), independent of 1 Although the 



below, implies that M = 0(M^ /2 ) 



C(e"^*), this will still be an 0(1) 
quantity with respect to v for all t < 0(1/ v). Thus, we may take T = 0(1/ v) in the statement of 
the Proposition. Gronwall's inequality then implies that 

(3.27) &(t) < <f> e (0)e~^\ 

We now use this estimate to control the norm of solutions. Throughout this estimate, C\,C2,K 
denote generic constants that are independent of any relevant quantity. Equation (3.25 ) and Remark 
|3.3| imply that 
(3.28) 



1 



\w\ 



\Wa 



Recall that 



l 2 4c 2 



Mix 



1 



v^l 3/2 



E 

1^0 



|C7^|| 2 < &(t) < \\wf+^C- 
l^l| 2 + V ^ll^f + 77^375 



v 'i - n2 3^ 

w x +-C 2 



1 



v^l 3/2 



\C l w\ 



Then, (3.28) and (3.27) imply 



and we obtain Theorem [TJ 



w 



(t)f x <Ke- M ^ l \ 



w 



2 

x ■ 



Proof of Proposition 3.4 We now prove the proposition via a series of lemmas, following closely 
the approach of Gallagher, Gallay and Nier. 

Define a partition of unity X+( x ) + X-( X ) = 1 where \± are C°° and non-negative with supp(x+) C 
[— 6, 7T+5], for some 5 small and positive, &ndx+(x) = 1 for all x 6 [8, w — 5], supp(x-) C [—tt — 5,5], 
and X-( x ) = 1 f° r & H x 6 [ — ^ + ^ ~~ £]■ F° r an Y u £ -^ 2 > define u± = x±^- 

Lemma 3.6. There exists a constant C > suc/i that for any u £ H , one has 

\\d x u\\ 2 > ^(Ru + || 2 + ||cV-|| 2 ) - C||n|| 2 



l^n+H 2 + \\d x u. 



3 
4 1 



> -\\d x u\\ 2 



C\\u\\ 
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Proof. First note that 

(3.29) J \d x u\ 2 = J d x u(x\ + X-)d x u = J d x ux+d x u + j d x u\ 2 -d x u. 
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One can write 



d x ux\d x u 



We can now bound 



(x+d x u)(x+d x u) 
(d x u^)(d x u+) - ( / {(d x x+)u)(x+d x u) + I (d x u^)((d x x+)u) 



(d x u + )((d x x+)u)\ < Ci\\u\\\\d x u\\ < -||c? x u|| 2 + C 2 ||ti|| 2 . 



Applying the same sort of estimates to the last term on the right hand side of (3.29), we obtain 

J \d x u\ 2 > {\\d x u + \\ 2 + ||^n_|| 2 ) - \\d x u\\ 2 - 4C 2 
from which the first inequality in the lemma follows. 



J 2 \\U | 



To prove the second inequality, we again use equation (3.29). Consider the first of the two integrals 
on the right hand side, which can be written 



{x+u x ){x+Ux) 



{d x u + - x'+u){d x u + - x'+u) 
J \d x u + \ 2 + J ( X ' + ) 2 \u\ 2 - 2ReJ (c\u^)x' + u 



< - j |d^+| 2 + c 3 / \u\ 2 , 



where the last inequality used the bounds J*(x+) 2 \u\ 2 < Ci||u|| 2 and |2Re f(d x ii+)x'+u\ < ^f\d x ii+\ 2 + 



C*2 / \u\ 2 - If we apply a similar estimate to the term J (x-U x )(x~u x ) in (3.29), we then obtain 
from which the second inequality in the Lemma follows. 



□ 



We now begin to assemble the pieces to prove Proposition |3.4| From the definition of C 
1 



-2ut 



■2 , ftii^.,ii2 ^ ^ no „.n2 , C 2 l foe" 2 



\d x u\\ z + —\\C l u\\ z > C 1 \\d x uf + 
2i> 



cosx u 



We now use Lemma 3.6 to bound this expression from below by 



2Ci... I|2 , C 2 l 2 (3 e 
-\\d x u + \\ +- 



2vt 



COSX lij 



W 2Cl ll* n2C 2 l 2 ^e 
+ I -^-\\dxu-\\ z + - 



-2vt 



COSX U- 



3 " v j \ 6 v 

Lemma 3.7. There exists M±, such that for any u+ with support in [—5, ir + 5], 



C 3 \\u\ 



-vt 



2Ci..- „ 2 C 2 / 2 /3 e-^ „ 2 J^J£jv%e- ||2 
— — 1 1 03.144. 1 1 "I ||cosxu + || > -= \\u+\\ 
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Proof. There exists D > such that for any —5 < x < tt + 5, cos 2 x > D{x — Tt/2) 2 . Thus 



2C 



1 no 
-\\o x u A 



+ 



C 2 l 2 (3 e 



-2ut 



COS X Ua 



> 



2d 



\d x u^. 



+ 



DC 2 l 2 p e 



-2vt 



-{x - tt/2) 2 \u + \ 



3 v J \ 3 v 

But, as pointed out in [61 Appendix A], this last integral is the quadratic form associated to the 
quantum mechanical harmonic oscillator. In more standard notation, 

H = -h 2 dl + u?x 2 =>• (Hu,u) L 2 (R) > hw(u,u) L 2 iR) . 

Transferring this to the current setting, we obtain the lemma. □ 



If we also apply Lemma 3.7 to u_, we obtain 



8 Iv Jv 



k+l| 2 + h-l| 2 )-C 3 | 



u 



Now apply Lemma 3.6 and find that 

2 , ^0 \, ,i|2 



-ut 



from which Proposition 3.4 follows if v is sufficiently small. ■ 
We conclude this section by briefly discussing how we might extend the theoretical results of this 



section to the full linearized operator (2.1). Villani's hypocoercivity method applies to operators 
of form 



(3.30) 



C = A*A + B 



where B is anti-symmetric. If we consider the linear operator (2.1) restricted to functions with 



fixed Fourier coefficient I, as we did for the approximation operator earlier in this section it can be 
written as 

(3.31) J~-£W£ = uAiW£ — iiae~ vt {smx)(l + A^" 1 )^ , 

where = (d 2 — £ 2 ). We can write vAp as A* A (as we did for the approximation linearization 
studied here), but the remainder of the operator is no longer anti-symmetric. 

However, note that on the space Ai defined in Section[2j (i.e. the complement of the "slow" modes), 
the operator (1 + A7 1 ) is a symmetric, strictly positive operator, and hence we can define its square 

root, \ 1 + A7 1 . Then, if we make the change of variables ue = \ 1 + Aj 1 W£, this transforms the 



operator L(W( to 
(3.32) 

The operator B e - 



Ciue = vAui — i£ae 



-ut 



1 + A7 1 sinWl + A7 1 



-iiae 



-vt 



l + A7 1 sinx A /l + A; 



is anti-symmetric, so we can attempt to 



apply Villani's theory to this new form of the operator. We note that in passing having written the 
operator in this form (i.e. the sum of a dissipative operator, plus an anti-symmetric operator), we 
immediately obtain that the bar states are stable because all eigenvalues of must have negative 
real part. However, this abstract argument gives only a weak estimate on how large these negative 
parts are - the hypocoercivity argument has the potential to give much stronger estimates on the 
location of these eigenvalues. 
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If we define 

\Jl + A7 1 cosx^/l + A7 1 

most of the preceding argument goes through without change. However the fact that [B e , C 1 -} 7^ 
causes problems at several points. We believe that this problem can be overcome by altering the 
functional 3> (t), and plan to address this point in future research. We note that similar ideas have 
recently been used by Deng [5] to study the linearization of the two-dimensional Navier-Stokes 
equation about the Oseen vortex solutions. 



(3.33) 



C e = [d x 



-iiae 



-vt 



4. Numerics 



The results of the previous section prove that for our approximation to the linearization about the 
bar states, solutions decay at a rate of at least 0(e~^' t ), much faster than the viscous decay rate 
of 0{e~ vt ). In this section we present numerical evidence that this rate of approach is actually 
optimal, not just for the approximate linearization studied in the previous section, but also for 



the full linearization (2.1). We numerically compute the eigenvalues of the linearization of the 
two-dimensional vorticity equation about the bar states, for a fixed value of the time t, and show 
that as we vary the viscosity parameter i/, the real part of these eigenvalues are negative and scale 
like C\fv, so long as we avoid the anomalous subspace. 

Recall that as we remarked earlier, if we expand the space L 2 ([— vr,7r] x [— tt, it]), with respect to 
the basis e lkx e lly , then the operator C{t) leaves invariant the subspaces with t fixed. Therefore we 
can study the eigenvalue asymptotics, by fixing i and studying the spectrum of the operator £i(t) 



defined in (3.31) for different values of £. We compute (approximately) the spectrum of Ct(t) by 
writing out its matrix representation with respect to {e tkx }^ = _ N and then using Mathematica to 
compute the eigenvalues of the resulting matrix. We have examined matrices with different values 
of N, up to N = 400 and we find consistent results for all values of iV provided N is large enough 
(bigger than about 30-40) to insure that we have good asymptotics. 

We present two examples of our numerical computations below, both with N = 400 and 1 = 2. 
Also, we fix the time t so that ae~ vt = 1. 

The first figure plots the logarithm of the absolute value of the real part of the least negative 
eigenvalue, as a function of the logarithm of u, for values of v = 0.005, 0.002, 0.001, 0.0005, 0.00025 
and 0.0001. Superimposed on these six data points is a line of slope 2, which corresponds to the 
real part of the eigenvalue scaling as \fv. The fit is quite good, especially for the smaller values of 
v. 



For our second piece of numerical evidence in support of the optimality of the yfu scaling of the 
attractivity of the bar states, we first compute the eigenvalues of linearization about the bar states, 
and then order the resulting eigenvalues according to decreasing real part. We consider the first 30 
eigenvalues in this ordering and plot the real part of each eigenvalue, divided by y/u. We repeat 
this procedure for three different values of u, u = 0.00025, v = 0.0001 and v = 0.00005. If the 
eigenvalues were of the form Ck y/v, j = k, . . . , 30, then these three data sets would fall on top of each 
other. As Figure [2] shows, this is nearly true, especially for the first few eigenvalues. As one moves 
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to eigenvalues further down the list, the smaller values of v overlap more than the larger values. 
This is to be expected - for larger values of k, the diagonal part of the linearization —u(l 2 + k 2 ) will 
dominate the behavior of the eigenvalue until v becomes quite small and the interaction between 
the diagonal and off-diagonal parts captured by the hypocoercivity method makes itself fully felt. 



5. Discussion 



We have analyzed the linearization of the two-dimensional vorticity equation (1.1) about the bar 
states. The bar states are "quasi-stationary" states since for small values of the viscosity, v, they 
decay on the very slow scale 0(e~ L>t ). We have proposed a new, dynamical systems based, expla- 
nation for the metastable behavior these states exhibit in numerical simulations. Our explanation 
begins by showing that a natural approximation to the linearization of the Navier-Stokes equations 
about the bar states exhibits decay on a much faster time scale than the viscous decay rate and 
also presenting numerical computations of the eigenvalues of the full linearization of the equation 
which are consistent with the decay rate we obtain rigorously for our approximation. 



Ideally, one would like to extend this linear result to the nonlinear equation ( 1.3 ). We are currently 
working on this extension, first of all by showing that our theoretical results on the approximate 
linearization can be extended to the full, linearized equations. The direct extension of the linear 
results to the full, nonlinear equation is complicated by the fact that the linearized operator is time 
dependent, so that we can't use ordinary spectral projections to project the nonlinear equation 
onto rapidly decaying and slowly decaying modes (the "anomalous" modes of the linear equation, 
described in Section [2]) as is done in the stable manifold theorem, for instance. 
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Real part of eigenvalue divided by square root of nu 
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Figure 2. The scaling behavior of the first thirty eigenvalues of the linearization 
about the bar states, with I = 2 and a 801 x 801 matrix approximation to the 
x-dependent part of the operator. In this figure, solid squares correspond to the 
eigenvalues with v = 0.00025, solid circles correspond to v = 0.0001 and clubs 
correspond to v = 0.00005. 



Another interesting direction to pursue is to analyze the stability of the dipole states, given explicitly 
by 

(5.1) UJ d,m (x, y, t) = e~ um2t [cos mx + cos my], £j d ' m (x,y,t) = e~ um2t [sin mx + sin my] 
with associated velocity fields 

u d ' m (x,y,t) = -e~» m2t ( ~ smm y ) , u d > m ( x ,y lt ) 
m \ sinmx / 

These states are also maximum entropy solutions of the Euler equations, and hence likely candidates 
for relevant quasi-stationary states for the 2D vorticity equation. In fact, these states seem to be 
more important than the bar states, at least for a square torus [3], as solutions with a larger class 
of initial data converge to them. See |18|, Figures 2a and 5]. However, we have not yet found an 
invariant subspace of solutions that converge rapidly to these states. 

If one linearizes about the dipole solution aw^ 1 , one finds 

(5.2) d t w = C d (t)w -v -Vw 

C d (t)w = vAw -ae- ut [sinxd y {l + A~ r ) - sinyd x (l + A^w. 



^ e -vmH I -cos my 
m \ cosmx , 
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In Fourier space, this operator becomes 



C t9 (t)w{k,l) 



-v{k 2 + l 2 )w(k,l) 



-ae 



-wt 



1 



(5.3) 



+ -ae 



-vt 



1 

' (A; + l) 2 + Z 2 
1 

k 2 + {l-l) 2 



(k 

w(k+ 1,1) 
w(k,l- 1) 



l) 2 + P 



w(k - 1, 1) 



k 2 + (l + l) 2 



,1 + 1) 



Although e x and e^ iy are still eigenfunctions, the adjoint eigenfunctions seem to be more difficult 
to determine. (An explicit one is cosx + cosy.) A similar change of variables can be used to 
transform the second part of the linear operator into a skew symmetric one: by defining w = 
+ A _1 u, the linear part of equation (5.2) becomes 

(5.4) u t 



i/Au + ae vt Qu, Q = \/l + A~ 1 \smyd x — smxd y ]\/l + A -1 , 



where Q is anti-symmetric. While we have some preliminary numerical data indicating that most 
modes again decay at a rate much faster than the viscous time scale, we have not yet proven a 
result similar to Theorem [T] for this operator. 
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